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ABSTRACT Drug-resistant extraintestinal pathogenic Escherichia coli (ExPEC) strains are the major cause of colisepticemia (coli- 
bacillosis), a condition that has become an increasing public health problem in recent years. ExPEC strains are characterized by 
high resistance to serum, which is otherwise highly toxic to most bacteria. To understand how these bacteria survive and grow in 
serum, we performed system-wide analyses of their response to serum, making a clear distinction between the responses to nu- 
tritional immunity and innate immunity. Thus, mild heat inactivation of serum destroys the immune complement and abolishes 
the bactericidal effect of serum (inactive serum), making it possible to examine nutritional immunity. We used a combination of 
deep RNA sequencing and proteomics in order to characterize ExPEC genes whose expression is affected by the nutritional stress 
of serum and by the immune complement. The major change in gene expression induced by serum — active and inactive — in- 
volved metabolic genes. In particular, the serum metabolic response is coordinated by three transcriptional regulators, Fur, 
BasR, and CysB. Fur alone was responsible for more than 80% of the serum-induced transcriptional response. Consistent with 
its role as a major serum response regulator, deletion of Fur renders the bacteria completely serum sensitive. These results high- 
light the role of metabolic adaptation in colisepticemia and virulence. 

IMPORTANCE Drug-resistant extraintestinal pathogenic Escherichia coli (ExPEC) strains have emerged as major pathogens, espe- 
cially in community- and hospital-acquired infections. These bacteria cause a large spectrum of syndromes, the most serious of 
which is septicemia, a condition with a high mortality rate. These bacterial strains are characterized by high resistance to serum, 
otherwise highly toxic to most bacteria. To understand the basis of this resistance, we carried out system-wide analyses of the 
response of ExPEC strains to serum by using proteomics and deep RNA sequencing. The major changes in gene expression in- 
duced by exposure to serum involved metabolic genes, not necessarily implicated in relation to virulence. One metabolic regula- 
tor — Fur — involved in iron metabolism was responsible for more than 80% of the serum-induced response, and its deletion ren- 
ders the bacteria completely serum sensitive. These results highlight the role of metabolic adaptation in virulence. 
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Septicemia, the presence of bacteria in the bloodstream, is 
among the leading causes of death worldwide, and its inci- 
dence is on the rise (1). In the past decade, there has been a rapid 
increase in the rates of hospitalization and mortality from severe 
sepsis, mainly because of the escalation of antibiotic resistance (2). 
While most bacteria are unable to endure the strong bactericidal 
effects of serum, several pathogens have evolved mechanisms that 
enable them to subvert the host defense systems and successfully 
survive in this hostile niche. 

In order to survive and even proliferate in serum, bacteria must 
overcome two major obstacles, the nutritional immunity and in- 
nate immunity of the host. Nutritional immunity is the process by 
which nutrients are kept in various storage molecules that make 
them unavailable to pathogens (3). Thus, an invading bacterium 
has to pass a metabolic barrier to survive. Iron sequestration by the 



host is the best-studied case of nutritional immunity. Although 
iron is an abundant nutrient in nature, serum contains very little 
free iron because iron is bound to storage molecules such as ferri- 
tin and hemosiderin. Therefore, it is not surprising that iron ac- 
quisition systems and receptors were found to play a pivotal role in 
the virulence of numerous pathogens. For instance, TonB, a pro- 
tein that provides energy for the transport of iron compounds, 
was found to be required for the pathogenicity of several Gram- 
negative bacteria (4-8). An additional factor that was shown to be 
important for growth in serum is the ability to synthesize nucleo- 
tides, as the inactivation of nucleotide biosynthesis genes was 
shown to hamper the growth of Escherichia coli, Salmonella en- 
terica serovar Typhimurium, and Bacillus anthracis in human se- 
rum (9). 

The second barrier that bacteria need to overcome to establish 
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FIG 1 ExPEC 078-9 is serum resistant. Cultures of E. coli strains 078-9 (A) and K-12 (B) were grown overnight and diluted in MOPS medium or in MOPS 
medium containing 40% serum, and growth curves were generated by measuring OD 600 with the BIO-TEK Eon platform. Viability was estimated by dropping 
10 fil of a culture treated with serum for 60 min on LB medium plates and incubating them overnight at 37°C (insets). S, serum; IS, inactive serum. 



sustained bacteremia is the one drawn by the immune system of 
the host. The complement system serves as the first line of defense 
against invading bacteria and acts on the outer membrane. In 
Gram-negative bacteria, the complement complex mediates di- 
rect killing by the formation of pores in the cell membrane. To 
avoid the highly bactericidal effect of serum, pathogens evolved 
structural features that inhibit complement-dependent killing. 
Many of these adaptations are in surface-exposed components, 
such as the outer membrane lipopolysaccharide (LPS) and the 
bacterial capsule (10). It was previously shown that variation in 
the length of the O antigen, as well as its type (E. coli for instance, 
has more than 190 types of O antigens), influences the level of 
serum resistance (11). Even a single nucleotide change that con- 
sequently results in a truncated gene coding for O antigen, is suf- 
ficient to turn a serum-resistant strain to a serum-sensitive one 
(12). However, not all of the adaptations are structural; some bac- 
teria can secrete complement-binding proteins and specific com- 
plement inhibitors. Staphylococcus aureus, for instance, employs 
an arsenal of proteins to interfere with and block the activation of 
the innate immune response ( 1 3 ) . Successful evasion of the innate 
immune system, coupled with proper metabolic adaptations, en- 
ables pathogenic bacteria to persist and multiply in the blood- 
stream, a process that often leads to sepsis and, if not treated suc- 
cessfully, death. 

Among the leading causative agents of sepsis are extraintestinal 
pathogenic E. coli (ExPEC) strains, which bring about high mor- 
bidity and mortality rates. E. coli infections also result in a heavy 
economic burden. ExPEC strains are characterized by the ability 
to resist the bactericidal effect of human serum, and this study was 
aimed at understanding the functional genomic basis of this resis- 
tance. To this end, we studied E. coli 078:H19 sequence type 88 
(ST88) isolate 789 (078-9), an ExPEC isolate able to grow, and 
even multiply, in serum. We performed comprehensive system- 
wide analyses of both its transcriptome and its proteome in re- 
sponse to serum. To distinguish between genes that respond to the 
metabolic challenge presented by serum and those that respond to 
the bactericidal effect of the complement system, we compared the 



effect of serum depleted of immune complement by heat treat- 
ment (inactive serum) with that of active serum. Here we show 
that the bacterial response to heat-inactivated serum included a 
change in the expression of a vast spectrum of known virulence 
determinants and iron acquisition genes, as well as transcriptional 
regulators. In contrast, the distinctive response to active serum 
(complement system) involved a much smaller group of genes, 
some of which are involved in the response to pH changes and 
others of which have unknown functions. One clear conclusion 
that emerged from these studies is that the key regulator of the 
response to serum — active and inactive — is the iron regulator Fur, 
which controls the expression of more than 80% of the serum- 
upregulated genes. 

RESULTS AND DISCUSSION 

ExPEC growth in serum. Septicemic bacteria are relatively resis- 
tant to the bactericidal effect of serum. We studied a septicemic 
strain, E. coli serotype 078:H19 ST88 isolate 789 (078-9), that is 
capable of growing in serum. The results presented in Fig. 1 A show 
the growth of E. coli 078-9 in 40% serum in comparison with that 
of a nonpathogenic E. coli strain (K-12) that undergoes lysis under 
the same conditions (Fig. IB). The serum component significant 
in growth inhibition is the immune complement, as its inactiva- 
tion by heat (inactive serum, 40 min at 56°C) enables the growth of 
E. coli K-12 (Fig. IB, inset). 

Transcriptome analyses of ExPEC genes upon exposure to 
serum. As pointed out, the response to serum is driven by two 
major stress conditions — nutritional stress and the bactericidal 
effect of complement. To study the nutritional stress, we investi- 
gated the response of ExPEC 078-9 to inactive serum. Because 
serum inactivation is achieved by mild heat conditions that de- 
stroy the complement system without affecting other properties 
such as pH and nutrient availability, the response to inactive se- 
rum represents the response to nutritional immunity. The exper- 
iments were carried out with cultures grown in minimal, salt- 
glucose, media. In these media, all of the biosynthetic genes are 
active and it is possible to study metabolic stress. 
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FIG 2 Tree map of differentially expressed genes in ExPEC 078-9 based on transcriptome analysis. Functional Voronoi tree maps of E. coli strain 078-9 were 
obtained from transcriptomes of bacteria exposed to serum for 15 min. Each cell represents one quantified transcript, and transcripts encoding functionally 
related proteins are subsumed in convex-shaped parental categories of increasing levels based on TIGR gene classification. The ratios of expression data were 
color coded by using a divergent color gradient. The bluish and orange tiles represent down- and upregulated pathways, respectively. Top panel: General 
representation of all functional pathways, levels 1 to 4 (left to right). Lower panel, left side: Transcription patterns in response to nutritional challenge. Lower 
panel, right side: Immune complement-induced gene fraction. 



To obtain transcriptional expression data for the complete ge- 
nome, we used the Illumina transcriptome sequencing (RNA-seq) 
platform. RNA was obtained from biological replicates of bacteria 
grown in minimal morpholinepropanesulfonic acid (MOPS) me- 
dium supplemented with glucose and with or without exposure to 
40% inactive serum. Total RNA was extracted and rRNA depleted 
(as described in Materials and Methods). Two independent bio- 
logical replicates and deep sequencing analyses were performed, 
and the data were analyzed by using strict criteria, as described in 
Materials and Methods. We examined only genes that were ex- 
pressed the same way in the two independent biological experi- 
ments by using a 2-fold cutoff. The expression data of selected 
genes were confirmed by quantitative reverse transcription-PCR. 

To obtain a comprehensive overview of the transcriptional re- 
sponse to serum, we compiled Voronoi tree maps (14, 15) from 
the quantitative transcriptomic data. In these tree maps, the tran- 
scripts are clustered according to their functions as obtained from 
the TIGR classification. Thus, functionally related elements are 
localized in close proximity to each other. The categories can then 



be subdivided sequentially down to the level of a single gene. A 
general representation of all functional pathways is shown in 
Fig. 2, upper panel. The log 2 ratios (i.e., treated/ control) of expres- 
sion data were color coded by using a divergent color gradient 
(from blue for downregulated to red for upregulated). The 
Voronoi tree maps of the transcriptome induced by inactive se- 
rum and the pathways derived from it are presented in Fig. 2, 
lower panel, left side, which shows that serum affects the expres- 
sion of a broad spectrum of genes with a variety of functions. Most 
of the serum-induced changes are metabolic; the upregulated 
pathways were mainly those of cations, iron-carrying compounds, 
amino acids, and amine biosynthesis, especially in the aspartate 
family. In contrast, pyruvate, glutamate, and aromatic amino acid 
biosynthesis pathways were downregulated. Other significantly 
downregulated genes included those for sulfur metabolism, the 
trichloroacetic acid cycle, and metabolism of polysaccharides. 
Most striking is the downregulation of biosynthetic genes involv- 
ing the synthesis of several amino acids, mainly cysteine. The 
downregulated genes (Table 1) appeared to be affected directly by 
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TABLE 1 Genes downregulated (at least 2-fold) in response to inactive serum 



Operon 


Gene(s) differentially expressed 


Function(s) 


Regulator(s) 


cysPUWAM 


cysP, cysU, cysW, cysA 


Sulfate/thiosulfate transporter 


CysB, H-NS 


L.ysLJL V Kj 


CySLJ , LySl \ . L ys ' - 


*siiIt;i+p -3 1\ pn v\ \!\Xva n cTf*v3 cp cmt\ ai\e*n\A\i\ — ci 1 lrot"P L"incicp 
OUlldlC dUCll Vly 111 dllfilCI dJ>C dllLl dUCllViyi JiLllldlC MlldSC 


CysB 


cysJIH 


CVJj , l-VJl-i LVJJi 


^ulntp ann ~\ -nnrwnnn9MpT*ivlvli;iilT:3tp i*pniir~tacpQ 

.JLlllllC dllU -t 1/11U5 IJllUdUCll VI V 13 LllldlC lt_ULll-ltljK_o 


CvsB THF 


cysK 


cysK 


Cysteine synthase A 


CysB 


cbl 


cbl 


CysB-like transcriptional activator 


CysB, NtrC 


amBCADTEF 


CLTtiB, cifttC, avnA, civttD, civnT, urnE 


Modification of LPS 


BasR 


yibD 


yibD 


Predicted glycosyltransferase 


BasR, NsrR, PhoB 


glnK-cwitB 




Nitrogen-regulatory protein 


Fur, NtrC, FNR, GadX 


ftnA 


ftnA 


Cytoplasmic ferritin iron storage protein 


Fur, H-NS 


basS-basR 


bass, basR 


Response to excess external iron 




n is 


ais 


Protein induced by aluminum 




yjdB 


yjdB 


Metal-dependent hydrolase 




ugd 


ugd 


UDP-glucose 6-dehydrogenase 




yeeED 


yeeE, yeeD 


Unknown 




yedEF 


yedF 


Unknown 




ydjN 


ydjN 


Unknown 




ydjO 


ydjO 


Unknown 





two major transcription factors, CysB and BasR. The cysB gene 
product downregulates 1 1 genes involved in cysteine metabolism 
in addition to that for the transcription factor Cbl, which is re- 
sponsible for sulfate starvation. The downregulated genes also in- 
clude the six members of the sugar nucleoside metabolism am 
operon (16-18) (arnB, arnC. arnA, arnD, arnT, arnE), which are 
under the control of the BasRS two-component system. The latter 
consists of two proteins that are involved in iron metabolism and 
protect against high concentrations of external iron (19, 20). Gen- 
erally, it appears that the downregulated genes encode mainly 
metabolic enzymes (aminotransferases, glycosyltransferases, and 
formyltransferases) . 

The metabolic changes in response to serum are characterized 
by a massive upregulation of operons attributed to various iron 
receptors, binding proteins, transporters, and siderophore genes 
(Table 2). We show that the vast majority of these genes are in- 
volved in the maintenance of iron equilibrium between the bac- 
terium and its surroundings. Interestingly, many of the upregu- 
lated genes have basic metabolic functions. For instance, the entire 
operons encoding iron hydroxamate (fhuACDB), iron dicitrate 
(fecABCDE), FeS cluster assembly protein (sufABCDSE), an elec- 
tron transport system (nrdHIEF), and enterobactin (fepA-entD) 
formed the scaffold of the ExPEC 078-9 iron response to serum. 
Notably, these operons are also found in commensal E. coli strains 
such as K-12. 

Most interesting is the upshift in the expression of bglj, a com- 
ponent of the yjjQ-bglJ transcription unit. These genes detoxify 
the a-oxoaldehyde methylglyoxal (MG), an intracellular metabo- 
lism product; overproduction of MG results in cell death. The 
gene products contain the helix-turn-helix motif typical of the 
LuxR transcript regulator family (21) and are therefore suspected 
of being transcription regulators. It was previously suggested that 
they;;Q gene of E. coli serotype 02:H5 (22) and S. enterica serovar 
Typhimurium (23) may be important for virulence in long-term 
systemic infections of mice. 

Another upregulated gene is hokE, a homologue of hok (host 
killing) that is part of the hok-sok toxin-antitoxin system, an asso- 
ciated gene member of the SOS regulon (24) . Five additional genes 
that are upregulated in the presence of serum have unknown func- 
tions (yrbL pqqL, 78901042, 78903310, and 78904620). 



Complement-dependent transcriptomic response. The genes 
responding to the bactericidal effect of serum could be identified 
by exposure to active serum and comparison of the transcriptome 
data to those obtained in the presence of inactive serum. This 
comparison revealed a small set of genes that were induced only in 
the presence of active serum — probably by the presence of com- 
plement. In contrast to the metabolic stress (Fig. 2, lower panel, 
left side), the complement-induced gene fraction is characterized 
by upregulation of the serine family, sulfur metabolism, and py- 
ruvate and glutamate biosynthesis (Fig. 2, lower panel, right side). 
The complement-induced response consists of genes involved in 
iron metabolism (Table 3) — the ferric enterobactin transport sys- 
tem gene fepE and the iron storage bacterioferritin gene bfr. The 
presence of complement also induces galP, a gene that codes for a 
transporter, a member of the major facilitator superfamily of 
transporters. The rest of the upregulated genes have unknown 
functions and code for the putative proteins YcfJ (25, 26), YpfG 
(27), and YgaC (28, 29). Exposure to complement repressed the 
transcription of the enzyme-encoding genes yfaO, sodB, sorC, and 
fimA in addition to one conserved hypothetical protein with high 
homology to iraM. The observation that FimA, the major E. coli 
fimbrial subunit protein, is downregulated in response to serum is 
in line with previous findings suggesting that under low-iron con- 
ditions, E. coli decreases the expression of type I fimbriae (30, 31). 
The iraM gene is induced under magnesium starvation and serves 
as a stabilizer of RpoS, the sigma factor of the general stress re- 
sponse (32, 33). The sodB gene, which encodes the Fe-dependent 
superoxide dismutase FeSOD, is also downregulated in the pres- 
ence of serum and is connected with iron metabolism (34). 

It should be noted that in addition to the genes whose expres- 
sion changes in response to complement, there are genes for other 
systems essential for complement survival, such as LPSs or cap- 
sule. These genes were studied in recent papers (10, 35) that 
looked at the effect of serum on an E. coli strain involved in urinary 
tract infections. Clearly, these surface components are essential 
for pathogenesis and survival in the bloodstream. Some of these 
genes are constitutive, and therefore, their expression is not af- 
fected by the presence of serum. However, the genes involved in 
the biosynthesis of the exopolysaccharide colanic acid are induced 
by serum (35) but were not detected in our system. These genes are 
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TABLE 2 Genes upregulated (at least 2-fold) in response to inactive serum 



Operon 


Gene(s) differentially expressed 


Function(s) 


Regulator(s) 


fhuACDB 


fhuA, fhuC, fhuD,fhuB 


Iron-hydroxamate transporter 


Fur 


fpftA-pntT) 


fpr>A pntD 

J C U 1\ , C 1 1 1 1 ' 


T* 1 n t p rn n a c 1 1 n 

l-jllLClWUtlCLlll 


Fur CRP 


fes-ybdZ-entF-fepE 


fes y ybdZ, entF 


Enterobactin 


Fur, CRP 


fepDGC 


fepG, fepC 


Enterobactin transporter 


Fur 


entCEBAH 


Cfl tC €11 tE CIltB 611 tA 6tl tH 


Enterobactin 


Fur, CRP 




sufA sufR sufC vifD wfS 


/"liiQtPi* mpmnlv Tti*r*tpin 

1 CO ClLlOlCl dJoClllUly L/L KJ LC111 


Fur THF TsrR NsrR ChrvR 

1 HI j 1111, loClv, 1 N jL lv, \_/A. y lv 


nrdHIEF 


ti t'dH h veil fi vdE ti vclF 


Electron transport system 


Fur, IscR, NrdR 


fecABCDE 


fecA, fecB, fecC, fecD, fecE 


Iron-dicitrate transporter 


Fur, CRP, PdhR, YjiE 


fecIR 


feci, fecR 


Signal transducer, sigma 19 




prhRD 


P\~nJi PVnJ~) 


l_jllClliy Ll CLlloVl L1C L±\J11 jyoLClll 


Fur 


fen A RC 


feoB 


Ferrous iron transporter 


Fur, NagC, FNR 


fhuE 


fhuE 


Ferric-rhodotorulic acid transporter 


Fur 


jhuE 


fhuF 


Ferric iron reductase 


Fur OyvR 

1 Lll , VyAVlV 


entS 


entS 


Enterobactin exporter 


Fur 


fetiB 


fepB 


Enterobactin transporter 


Fur, RutR 


ybiX 


ybiX 


Fe(II)-dependent oxygenase 


Fur 






Catecholate siderophore receptor 




bfd-bfr 


bfd 


Bacterio ferritin, iron storage, detoxification 


Fur" 






Unknown 


Fur a Mar A 


ybtA 


ybtA 


Yersiniabactin transcriptional regulator 


Fur n 


iron 


iroN 


Outer membrane siderophore receptor 


Fur rt 


vddAR 

V if 1 1/1 i) 


vddA vddR 


AAiilfiHviio trancnnrtpi* 

1V1 LllLldlLli^ LldlloL'UllCl 


1 Lll \ y LiiAl \ j 


efeU 


efeU 


Ppvitihq i nil tvatiQnni'tPT' 

1 Ll 1 U Llj 1LJ11 Ll til 1 . » |_ \ ) 1 LCI 


Piii* n f^nvT? 

1 Ll 1 , L" A.1X 


efeOB 




JTCIIULIS 1U11 LI d-HSUUl LCI 




yjjQ-bgl] 


hall 


DNA-binding transcriptional activator 


rl-INo, LcllU 


hokE 


hokE 


Toxic polypeptide 


LexA 


ybil 


ybil 


Unknown 




yrbL 


yrbL 


Unknown 


BasR, PhoP 


pqqL 


pqqL 


Predicted peptidase 




78900613 


78900613 


Ferric enterochelin esterase 


Fur" 


78901042 


78901042 


Phage-like element 




78903310 


78903310 


Unknown 




78904620 


78904620 


Unknown 





,! Fur binding box computationally detected in promoter region. 



under the control of the Res regulator, whose level is greatly af- 
fected by the growth medium (36). Thus, the experiments of Mi- 
ajlovic et al. were conducted with cultures grown in Luria-Bertani 
(LB) medium, where the expression levels of the Res-regulated 
genes are low. Our experiments were carried out with cultures 
grown in minimal medium, where the Res-regulated genes are 
already induced and the addition of serum probably does not 
result in significant further induction. 



Proteome analyses of the serum response. The effect of serum 
on gene expression was determined at the proteome level by using 
mass spectrometry (MS) -based proteomics as described in Mate- 
rials and Methods. As in the transcriptome experiments, the pro- 
teins were extracted from cultures exposed to active or inactive 
serum and an untreated reference culture. A total of 1 ,037 proteins 
were identified (see Appendix SA in the supplemental material), 
and many of them could be detected only under one condi- 



TABLE 3 Complement-dependent genes 
Genes differentially 



Operon 


expressed 


Function(s) 


Response to serum'' 


Regulation 


fes-ybdZ-entF-fepE 


fepE 


Enterobactin 


Up 


Fur, CRP 


ygaC 


ygaC 


Response to cytoplasmic pH stress 


Up 


Fur 


ypfC 


ypfG 


Unknown 


Up 


Fur" 


ycfJ 


y<$ 


Biofilm formation 


Up 




bfd-bfr 


bfr 


Bacterioferritin, iron storage, detoxification 


Up 




galP 


galP 


Galactose-proton symport of transport system 


Up 


CRP, GalR, GalS, NagC 


78901228 


78901228 


Enhancement of lycopene biosynthesis protein 1 


Down 


ArcA, FNR, PhoP 


fimAICDFGH 


fimA 


Major type 1 fimbrial subunit 


Down 


H-NS, IHF, Lrp 


sodB 


sodB 


Superoxide dismutase 


Down 


NsrR, CRP, IHF 


yfaO 


yfaO 


Putative nudix hydrolase YfaO 


Down 




sorC 


sorC 


Transcriptional regulator of sorbitol uptake and 
utilization 


Down 





a Fur binding box computationally detected in promoter region. 
b Up, upregulation. Down, downregulation. 
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FIG 3 Tree map of differentially expressed genes in ExPEC 078-9 based on proteome analysis. The tree map shown is as described in the legend to Fig 
is based on proteomic data. 



2 but 



tion — in the presence of serum or in the untreated culture. Two 
hundred proteins were significantly differentially expressed, at 
least 2-fold, in the presence of serum (113 were upregulated, and 
87 were downregulated). Of these proteins, 155 were also differ- 
entially expressed in response to inactive serum (70 were upregu- 
lated, and 85 were downregulated). The proteins that are present 
at very low levels in the presence of serum are probably encoded by 
genes that were significantly repressed in the presence of serum. 
However, it is possible that the level of some of them goes down in 
the presence of serum because of degradation. In an effort to 
achieve a global view of the functions altered in response to innate 
or nutritional immunity, we mapped all of the proteins to their 
corresponding pathways on the basis of the Voronoi tree map 
representation (Fig. 3, top panel). 

The group of proteins that were induced by serum included a 
number of conserved hypothetical proteins. However, many of 
the proteins that were present in very low concentrations without 
serum and induced by serum (inactive and active) were involved 
in iron acquisition. This finding is not surprising, as serum is 
devoid of free iron because of the presence of human iron-binding 
proteins. There were also proteins that could not be detected once 
serum (active or inactive) was added. This group includes proteins 
that participate in various biosynthesis pathways, the am operon 



{arnC, arnA, and arnD [16-18]), transcription regulatory factor 
BasR (19, 20), and conserved hypothetical proteins. 

A comparison of the genes differentially expressed in the pres- 
ence of active serum versus inactive serum is shown in the bottom 
part of Fig. 3, right side. This group contained only a few proteins, 
as most of the responding proteins behave the same in active and 
inactive serum. The unique group of complement-induced pro- 
teins contained proteins involved in energy metabolism, amino 
acid synthesis, and the biosynthesis of cofactors, purines pyrimi- 
dine, and nucleosides. This group also contained proteins in- 
volved in cell envelope maintenance (MurD, MurE, M reB, GlmU, 
Air, GalU). The latter observation is in line with the knowledge 
that the innate immune system causes membrane stress and that 
the integrity of the bacterial cell envelope is crucial for survival in 
serum and circumvention of the host's membrane attack complex 
(37-39). 

For the complete list of proteins whose expression depends on 
the presence of serum, see Appendix SA in the supplemental ma- 
terial. 

Comparison of mRNA levels and protein abundance. An 

analysis of the mRNA levels and protein concentrations in re- 
sponse to exposure to inactive and activated serum (Fig. 4A and B, 
respectively) shows that under both conditions, there is a linear 
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FIG 4 Comparative analysis of RNA and protein profiles of E. coli 078-9 in 
response to human serum. The results show the correlation between changes 
in protein (log, fmol/ng ratio) and transcript (log, RPKM ratio) levels in 
bacteria grown in 40% serum over bacteria grown in MOPS minimal medium. 
A 2-fold change between the conditions was used as the cutoff for up- and 
downregulation of RNA and protein levels. (A) After exposure to inactive 
serum. (B) After exposure to active serum. 



relationship between RNA levels and protein abundance (active 
serum, R 2 = 0.3 and Spearman's rank correlation coefficient = 
0.45; inactivated serum, R 2 = 0.25 and Spearman's rank correla- 
tion coefficient = 0.47). These results are consistent with previous 
observations that mRNA levels explain one- to two-thirds of the 
variance in protein levels (40). 

Fur is the major regulator of the serum response. In an at- 
tempt to discover the regulatory network underlying serum adap- 
tation, we systematically analyzed which transcription factor con- 
trols each of the differentially expressed genes involved in ExPEC's 
serum response. To this end, we used a literature survey and Regu- 
lonDB, a database of the regulatory interactions of E. coli K-12 
(41). As shown in Table 2, in accordance with the proteome re- 
sults, Fur appears to be the key regulator of serum adaptation, 



g i.o- 



A 



C 



Position 

sequence 

...TTACGCTATTGCAAATGCAAATAGTTATCAATAATATTA.. 
...TTTTCTTTGCCGTAATGATAATCGATATCACTGCGATTT.. 
...TTTTCTTTTGCTAAATGATATTAATTATCATTCTTATAT.. 
...CATTTGAGAAAATAATGATTACCATTCCCATTTATAACA.. 
...GTAAATCATTGCCATAGATAATCATTGGATGGCACAATC. 
...GAGAATTAATGTGAATAATAACCATTATCAATAACTCAT.. 
...TAATGGACTGAAAAATAACAATAATATTTCGTTGTTTAT.. 
...TTTCCGTAATTGAAATAAGAACTATTTTCATTTATTTTT.. 



Gene 

78900613 
efeU 
iroN 
yncE 
yddA 
ybtA 
ypfG 
bfd 

FIG 5 The Fur binding box motif and predicted E. coli 078-9 genes with Fur 
binding boxes in their promoter regions. The Fur binding box 19-bp consen- 
sus sequence was created in WebLogo (upper part). Eight 078-9 gene se- 
quences containing a Fur binding box (highlighted) that were computationally 
derived from the PRODORIC database are shown at the bottom. 



controlling more than 70% of all upregulated genes. Fur, the ferric 
uptake regulator, serves as the main iron regulator in a wide range 
of bacterial species (29, 42-51). Fur was found to play a role in the 
virulence of many bacteria, including Staphylococcus aureus, 
Vibrio cholerae, Neisseria species, Yersinia species, and E. coli (52, 
53). Fur repression occurs under iron-depleted conditions and is 
achieved by binding of Fur dimers to Fur boxes of target promot- 
ers. In the absence of free iron, the repression is removed and the 
genes are transcribed, as occurs in serum, an iron-limiting envi- 
ronment. 

While RegulonDB provides comprehensive data on transcrip- 
tion factor-gene interactions in E. coli K-12, it cannot detect new 
regulatory interactions that appear in other strains. To overcome 
this obstacle, we scanned the region upstream of each of the 
serum-responsive genes to identify the presence of the 19-bp Fur 
binding consensus box (Fig. 5, upper part) by using weight matri- 
ces taken from PRODORIC (54). The Fur binding box is com- 
posed of three adjacent hexamers that bind dimeric Fur (55). Fur 
binding boxes were detected in seven additional genes that were 
upregulated in response to inactive serum, efeU, iroN, yncE, ybtA, 
bfd, yddA, and the ferric enterochelin esterase-encoding gene 
78900613 (Fig. 5, lower part), and in one gene, ypfG, that was 
induced only in the presence of complement. As expected, most of 
these genes are attributed to iron transporters, acquisition sys- 
tems, and storage. It may be significant that two out of the eight 
genes found to have a Fur box, iroN and ybtA, were acquired by 
lateral gene transfer from other pathogens (S. enterica and Yersinia 
pestis, respectively) (56, 57). Thus, all in all, Fur is predicted to 
control more than 80% of the genes involved in the adaptation of 
E. coli to serum (Fig. 6 and Table 2). 

Fur is essential for survival in serum. As Fur is clearly a major 
regulator of growth in serum, we constructed a deletion of the fur 
gene to study its effect on survival in serum. The results presented 
in Fig. 7 indicate that the removal of Fur had no effect on growth 
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study 

FIG 6 The majority of the genes upregulated in response to serum are controlled by Fur. Genes upregulated following exposure to inactive serum that contain 
a Fur binding box (blue) and genes upregulated under another transcription factor (red) are shown. Asterisks indicate genes predicted by this study to have a Fur 
binding box. 



in minimal medium. In contrast, the deletion-containing bacteria 
grew very poorly in the presence of serum. Growth in serum was 
restored by complementation with the fur gene on a plasmid. Al- 
though there was an induction of many proteins involved in iron 
metabolism, the concentration of Fur was unchanged upon expo- 
sure to serum. These results further support the notion that the 
effects of Fur are not solely concentration dependent and there 



must be additional factors involved, such as iron-dependent 
changes in its activity (58, 59). 

Our results reflect the complex role of Fur in bacterial metab- 
olism. As Fur is the regulator of iron metabolism, its removal 
would disrupt the delicate balance of the function involved in this 
important system. Actually, one would expect that a fur deletion 
would have a significant effect under all growth conditions. Yet, 
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FIG 7 Fur is essential for growth of E. coli 078-9 in serum. Bacterial cultures were grown as described in the legend to Fig. 1. Strain 078-9 A/ur grew in MOPS 
medium (continuous line) or in the presence of 40% serum (dotted line). Growth arrest due to serum exposure was restored in strain 078-9 A/ur/pBAD/ur 
(dashed line). Viability was estimated as described in the legend to Fig. 1, and the results are shown on the right. S, serum; IS, inactive serum. 
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TABLE 4 Plasmids and strains used in this study 



Plasmid or strain Description Reference 
Plasmids 

pKD4 Template for kanamycin resistance cassette 76 

pBAD24 Amp 1 , arabinose-inducible plasmid 69 

pBAD/ur Amp 1 , arabinose-inducible plasmid containing fur ORF This study 

E. coli strains 

MG1655 Wild-type K-12 strain 

078-9 Wild-type ExPEC 078 strain isolated from a turkey with sepsis 

078-9 Afur This study 

078-9 Afur/pBADfur This study 



the bacteria appear to be able to control growth in the absence of 
Fur but are not capable of doing so under extreme iron concen- 
tration changes, such as in the presence of serum. The effect of fur 
deletion is not restricted to E. coli, as it was also previously shown 
in Neisseria meningitidis, deletion mutants of which acquired se- 
rum sensitivity (60). 

Conclusions. Here we present data, obtained from proteomic 
and transcriptomic experiments, on the global effect of exposure 
to serum on septicemic E. coli. We focused on genes required to 
overcome the nutritional immunity of serum and therefore con- 
ducted our studies with minimal media, in which the metabolic 
genes are active. Under such conditions, it is possible to see the 
serum-induced metabolic changes and the regulatory networks 
that control them. The most striking finding is the identification 
of Fur as the major regulator during serum exposure. Thus, >80% 
of the upregulated genes are under the control of Fur, the major 
regulator of iron metabolism (61, 62). 

Genes that encode virulence factors are essential for pathogen- 
esis. However, recent studies indicate that bacterial virulence is 
highly dependent on the regulation of metabolic genes (63-66), 
which are no less important than genes that encode virulence fac- 
tors. The results presented here strongly support this notion and 
indicate that the major changes in bacterial gene expression upon 
exposure to serum involve metabolic genes, and these changes are 
essential for survival in serum. 

MATERIALS AND METHODS 

Bacterial strains and growth media. All of the E. coli strain and plasmids 
used in this study are listed in Table 4. Unless stated otherwise, all E. coli 
strains were grown while shaking at 30°C in minimal defined MOPS me- 
dium (with 0.2% glucose added). When required, antibiotics were added 
to the medium (ampicillin at 300 u.g/ml, kanamycin at 50 u.g/ml, chlor- 
amphenicol at 20 u.g/ml). To determine growth curves, log-phase cultures 
were diluted to an optical density at 600 nm (OD 6U0 ) of 0.04 as determined 
with a BioTek Eon plate reader, and turbidity at 600 nm was measured 
every 20 min. Viable counts were determined after 1 h of incubation with 
serum or inactive serum by plating 10-u.l drops of increasing dilutions on 
LB medium plates and incubating them overnight. 

RNA isolation and purification. Overnight cultures were diluted to 
an OD 600 of 0.04 and allowed to grow to an OD 600 of 0.3; cells were then 
introduced to serum, inactive serum, or saline to a final volume of 40%. 
Cells were pelleted by 10 min of centrifugation at 3,000 rpm at 4°C. Pellets 
were resuspended in 6 ui of freshly prepared lysozyme solution ( 1 mg/ml) 
in 10 mM Tris-HCl, pH 8.0, and incubated at 37°C for 7 min with occa- 
sional mixing. One milliliter of TriReagent was added to each sample, and 
the homogenate was stored at room temperature for 5 min. A 200-u.l 
volume of chloroform was added for 15 s of incubation, and the resulting 
mixture was stored at room temperature for 2 to 15 min and centrifuged 
at 12,000 rpm for 15 min at 4°C. RNA was precipitated from the aqueous 



phase by mixing with isopropanol (0.5 ml/ml of the TriReagent used for 
the initial homogenization), stored at room temperature for 5 to 10 min, 
and centrifuged at 12,000 rpm for 8 min at 4 to 25°C. The RNA pellet was 
washed with 75% ethanol, centrifuged at 7,500 rpm for 5 min, and air 
dried. The pellet was dissolved in 30 to 50 u.1 of ultrapure water (or treated 
with diethyl pyrocarbonate). rRNA was depleted by using the MICROB- 
Express kit (Ambion) according to the manufacturer's protocol. DNase 
treatment was done according to the manufacturer protocol with the Am- 
bion TURBO DNA-free kit. 

RNA-seq. Removal of 16S and 23S rRNAs from total RNA was per- 
formed with the MICROBExpress bacterial mRNA purification kit (Am- 
bion) according to manufacturer's protocol (67). RNA concentrations 
and quality were determined by using the Bioanalyzer 2100 (Agilent) . For 
whole-transcriptome sequencing (RNA-seq), cDNA libraries were pre- 
pared by the Illumina mRNA-seq TruSeq protocol according to the man- 
ufacturer's protocol without the poly(A) isolation stage. Sequencing was 
performed with an Illumina HiSeq 2500 sequencing machine. Sequencing 
generated 20 to 28 million raw reads per sample, of which 2 to 6 million 
were mapped to protein-encoding genes. The RNA-seq data were ana- 
lyzed by a standard protocol based on the number of reads per kilobase of 
coding sequence length per million reads analyzed (RPKM) (68). 

Construction of deletions and recombinant plasmids. All site- 
specific gene knockouts were obtained as described by Datsenko and 
Wanner (69). Briefly, competent wild-type 078-9 bacteria were trans- 
formed with plasmid pKD46. The transformants were grown in 
ampicillin-containing LB medium, induced with arabinose, made com- 
petent for electroporation, and stored at — 70°C until used. A linear PCR 
product was made on the template of a kanamycin resistance cassette 
flanked by FLP recognition target (FRT) sequences from the pKD4 plas- 
mid according to the region to be deleted. The primers were designed to 
contain 36 nucleotides from the flanking region of the sequence to be 
deleted from strain 078-9. Kanamycin-resistant recombinants were 
screened by means of colony PCR with primers k2 and kt. Unless stated 
otherwise, the pKD46 plasmid was cured by growth on LB medium at 
42°C. The resulting bacteria were transformed with pCP20 and grown on 
LB medium at 42°C to promote both FRT-specific recombination and 
curing of the plasmid. The final deletion was verified by sequencing. Re- 
combinant plasmids for complementation were constructed by cloning 
PCR-amplified DNA fragments into a pBAD24 vector with Fast-Link 
DNA ligase (Epicentre). The ligation products were transformed into 
competent 078-9 bacteria, and the resulting clones were screened by col- 
ony PCR. 

Human serum. Sterile and filtered human male AB plasma was ac- 
quired from Sigma-Aldrich. Human serum was inactivated by incubation 
for 40 min at 56°C to eliminate all of the active immune complement 
system. 

Clustering and enrichment analysis. E. coli 078-9 is still available in a 
roughly annotated version. HMMER (70) was used for the prediction of 
gene functions according the TIGRFAM classification. For this purpose, 
HMMER used the Hidden Markov file containing the definitions derived 
from sequences of already assigned TIGR protein families from other 
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bacteria, which can be downloaded from the website ftp://ftp.jcvi.org/pub 
/data/TIGRFAMs. In the resulting list, locus tags (RACOXXXX) were 
assigned to one or more TIGRFAMs, which were assigned to (metabolic) 
subroles, which were then summarized in hierarchically parent main 
roles. TIGRFAMs that were assigned to generic categories such as hypo- 
thetical proteins, general, and so on were, if possible, manually reclassified 
into more-specific categories. By using Voronoi tree maps (14, 15), all 
functionally assigned genes were shown in a hierarchically organized and 
space-filling manner, as shown in Fig. 2. The space-filling approach di- 
vides the two-dimensional plane into subareas according the main rules, 
which are then subdivided into subsubareas of the subrules and so forth. 
The polygons on the last level represent all of the functionally assigned 
genes. Transcriptional data were mapped by using a divergent color gra- 
dient starting with gray in the middle and ending with blue (for repres- 
sion) and orange (for induction) on the sides. 

Identification of Fur binding sites. Fur binding site position weight 
matrices (PWMs) were downloaded from PRODORIC version 8.9 (54). 
The score of each sequence is the log of the ratio of the likelihood of the 
sequence given the PWM model to the likelihood of the sequence given a 
background model. The background model assumes that at each position, 
the probability of each character equals its frequency in a concatenated 
regulatory sequence of 078-9. A significant hit is defined as a score of >4.6. 
The position-specific scoring matrix for the consensus logo was created 
with WebLogo. 

Gel-free absolute proteome quantification. E. coli 078-9 was grown 
in 50 ml of MOPS medium in shaking flasks under agitation at 37°C. At an 
OD 500 of 0.5, human serum, inactivated human serum, or saline was 
added to a final concentration of 40% and the bacteria were incubated for 
45 min. Cells were harvested by centrifugation (8,000 X g, 15 min, 4°C). 
The resulting cell pellet was resuspended and washed twice in 50 mM 
triethylammonium bicarbonate. Cell disruption was carried out in a Ri- 
boLyzer (Thermo Fisher Scientific, Waltham, MA) for 2 X 30 s at 
6,800 rpm. Two centrifugation steps were used to remove glass beads and 
cell debris. The supernatant was used for the following analysis. 

Protein extracts (250 u.g) were digested with trypsin as described by 
Muntel et al. (71); this was followed by a desalting step via stage tips using 
a standard protocol described by Rappsilber et al. (72). For absolute quan- 
tification, an internal standard protein (tryptic digest of alcohol dehydro- 
genase; Waters, Milford, MA) was used to spike the samples at a final 
concentration of 50 fmol/ fx\. 

Data acquisition via LC-IMS E setup. Peptide samples were analyzed 
with a nanoACQUITY ultraperformance liquid chromatography (UPLC) 
system (Waters) coupled to a Synapt G2 mass spectrometer (Waters). For 
each sample condition, three biological replicates were acquired with 
three technical replicates each. The samples were loaded directly onto the 
analytical column (nanoACQUITY UPLC column, BEH130 C 1S , 1.7 ixm, 
75 by 200 mm; Waters) with 99% buffer A (1% acetic acid), 1% buffer B 
(1% acetic acid in acetonitrile) at a flow rate of 300 nl/min within 35 min. 
Separation of peptides was achieved in 165 min by applying the following 
gradient: to 18% buffer B in 102 min, to 26% buffer B in 22 min, to 99% 
buffer B in 16 min, at 99% buffer B for 10 min, and equilibration for 
15 min with 99% buffer A. The Synapt G2 mass spectrometer (Waters) 
was equipped with a NanoLockSpray source in positive mode and oper- 
ated with the MassLynx V 4.1 software (Waters). The analyzer was set to 
resolution mode and operated with the ion mobility separation. For anal- 
ysis of peptide ions, a mass range of 50 to 2,000 Da Th was used, low- 
energy (MS) scans were set to 4 V (collision energy [CE]), for elevated- 
energy ion MS (IMS E ) scans, the CE was ramped from 25 to 40 V, the scan 
time was set to 1 s, the wave velocity was ramped from 1 ,000 to 400 m/s, 
and the wave height was set to 40 V. For lock mass correction, [Glul]- 
fibrinopeptide B solution (m/z 785.8426 Da, 500 fmol/iil in 50% [vol/vol] 
acetonitrile- 0.1% [vol/vol] formic acid) was constantly infused at a flow 
rate of 500 nl/min and scans were acquired at intervals of 30 s. 

Analysis of LC-IMS E data. For identification and quantification of 
proteins, raw data were imported into ProteinLynx Global Server (PLGS) 



2.5.3 and processed via an Apex3D algorithm. The following processing 
parameters were used. The chromatographic peak width and MS time-of- 
flight resolution were set to automatic, the lock mass for charge 2 was 
785.8426 Da/e, the lock mass window was 0.25 Da, the low-energy thresh- 
old was 250 counts, the elevated-energy threshold was 30 counts, the 
retention time window was set to automatic, and the precursor/fragment 
ion cluster intensity threshold was 1,000 counts. A database search was 
carried out by the ion accounting algorithm implemented in PLGS (67) by 
using a database that consisted of 4,628 E. coli 078-9 entries in a raw 
database (73-75; Huja et al., unpublished data) and 51,767 human entries 
in a human reviewed database uploaded from UniProtKB. This database 
was complemented with common laboratory contaminants and the yeast 
ADH1 sequence. The following parameters were used for positive protein 
identification. The peptide tolerance was set to automatic; the protein 
tolerance was set to automatic; the minimum number of fragment ion 
matches per peptide was set to 1 ; the minimum number of fragment ion 
matches per protein was set to 5; the minimum number of peptide 
matches per protein was set to 1 ; the primary digest reagent was trypsin; 
the number of missed cleavages was 2; the variable modifications were 
carbamidomethylation C ( + 57.0215), deamidation N, Q (+0.9840), and 
oxidation M ( + 15.9949); the false-discovery rate (FDR) was 5%; and the 
calibration protein was yeast ADH1. 

All of the proteins identified in only one or two of the nine replicates 
(three technical replicates of each of the three biological replicates) of a 
sample were discarded. This filter procedure has revealed protein results 
with an FDR of <2% on the protein level. Proteins were considered to 
have significant expression differences if they showed at least a 2-fold 
change in serum or inactivated serum relative to the MOPS control. Only 
genes that were detected by both LC-IMS E and RNA-seq were used for the 
comparison of mRNA levels and protein abundance. 

SUPPLEMENTAL MATERIAL 

Supplemental material for this article may be found at http://mbio.asm.org 
/lookup/suppl/doi:10.1128/mBio.01460-14/-/DCSupplemental. 
Appendix SA, XLSX file, 0.1 MB. 
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